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We present a numerical investigation of shot noise suppression in mesoscopic cavities and an 
intuitive semiclassical explanation of the behavior observed in the presence of an orthogonal magnetic 
held. In particular, we conclude that the decrease of shot noise for increasing magnetic held is the 
result of the interplay between the diameter of classical cyclotron orbits and the width of the 
apertures defining the cavity. Good agreement with published experimental results is obtained, 
without the need of introducing htting parameters. 
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In the recent literature, the topic of shot noise sup- 
pression in mesoscopic structures has received significant 
attention, as a result of the formulation of several theo- 
retical predictions glbbbb and the subsequent exper- 
imental confirmation !g J7] of some of such predictions. 
All of the shot noise suppression phenomena are the re- 
sult of correlations between charge carriers that decrease 
the variance of the random process corresponding to the 
elementary charges crossing the device: such correlations 
may result either from Fermi exclusion or from Coulomb 
interaction, and for their investigation powerful theoreti- 
cal methods have been developed, ranging from Random 
Matrix Theory fg] to semiclassical techniques [q]. The 
convergence, in terms of shot noise suppression, between 
the results of quantum and semiclassical approaches has 
been explained, using a voltage probe technique, by van 
Langen and Biittiker [Till , who have shown that dephas- 
ing phenomena do not have effects on the noise power 
spectral density. This result appears to be valid also in 
the presence of non DC bias [Of • A detailed discussion of 
the fundamentals of shot noise in mesoscopic conductors 
can be found in ref. 0] ■ 

Particular interest has been raised by the so-called 
“chaotic cavities” , mesoscopic regions delimited by in- 
put and output constrictions that are much smaller than 
the cavities themselves. Jalabert et al. |2j showed that in 
the case of symmetric apertures noise is suppressed down 
to 1/4 of its full shot value: this theory received experi- 
mental confirmation Jjj in 2001, and, more recently, the 
noise behavior of a chaotic cavity in an orthogonal mag- 
netic field has been measured [l3j , observing a somewhat 
linear reduction of the Fano factor as the magnetic field is 
increased. In ref. the authors motivate the decrease 
of the Fano factor as the magnetic field is increased with 
the reduction of the portion of the cavity area explored 
by the electrons, and rely on a fitting parameter, the 
quantum scattering time, which is also used to explain 
the behavior of the Fano factor with no magnetic field as 
the apertures are made wider. They also attribute the 
differences between the variation of the conductance and 
that of the Fano factor to a transition between “quan- 
tum” and “classical” chaos. 

We provide an alternative, intuitive interpretation of 



the behavior in the presence of a magnetic field, based on 
the comparison between the classical cyclotron diameter 
and the constriction width, instead of the cavity dimen- 
sions, which is supported by the results of the numerical 
simulations that we present in the following. Further- 
more, we argue that there is no significant contribution 
of “classical” chaos to the observed phenomena, in par- 
ticular we suggest a different explanation of the preserved 
additivity of the constriction resistances up to 40 propa- 
gating modes, in the case of no magnetic field. Such an 
explanation is based on assuming some degree of nonide- 
ality in the experiment, specifically partial thermaliza- 
tion to the lattice temperature of the electron gas in the 
cavity. 

In most of our calculations, we consider model cavi- 
ties defined by hard walls and with a rectangular shape 
(see inset of fig. [T[, in order to keep the computational 
time within reasonable limits, running a few checks for 
structures with different geometries: only in few partic- 
ular cases small differences were observed, which will be 
detailed in the discussion of the results. We wish to 
point out that a classically chaotic shape is not needed 
to achieve the known results for shot noise suppression. 
This conclusion cannot be found in unequivocal and ex- 
plicit form in the existing literature, except for the very 
recent paper by Aigner et al. M , who have shown that 
the distribution of transmission eigenvalues (and there- 
fore the Fano factor) in a cavity with narrow constrictions 
does not depend appreciably on its shape. 

However, the theoretical basis needed to reach such 
a conclusion had already been developed by several au- 
thors. For example, we can apply to a single cavity, in 
the absence of a magnetic field, the analytic procedure 
discussed by Oberholzer et al. |T(| for a series of cavities. 
The cavity is supposed to act as an elastic quasi- reservoir, 
effectively decoupling the two constrictions: by balanc- 
ing the incoming and outgoing electron fluxes at each 
energy [^, one obtains, in the hypothesis of symmetric 
and narrow constrictions, the well-known result of 1/4 
for the Fano factor. 

We wish to point out that to arrive at this result, no 
hypothesis of a classically chaotic shape is needed: the 
discontinuities represented by the transitions between the 
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FIG. 1: Ratio of the cavity resistance to the sum of the con- 
striction resistances vs. the number of modes propagating in 
the constrictions (upper panel): the solid line is for our rect- 
angular cavity and the dotted line for a stadium-shaped cavity 
with a central constant-width region 1 /im long. Numerical 
(empty circles) and experimental (empty squares, from Ober- 
holzer et al. Hfl l results for the Fano factor as a function 
of the number of modes propagating in the constrictions, for 
no applied magnetic field (lower panel). The inset contains a 
sketch of the model cavity. 



cavity and the leads produce the mode mixing needed to 
generate a uniform occupancy of all the states at the 
same energy ( i.e . the property defining an elastic quasi- 
reservoir), which, along with symmetry and integer trans- 
mission of the constrictions, represents the essential con- 
dition needed to achieve the 1/4 suppression factor jlH . 

This is the reason why we prefer to define the sub- 
ject of our investigation as “mesoscopic cavities,” rather 
than “chaotic cavities,” in order to avoid possible confu- 
sion with classical chaotic dynamics. Quantum diffrac- 
tion occurring only at the input and output openings 
leads to multiple trajectories that create the equivalent 
of a chaotic behavior. 

Our numerical technique is based on the evaluation of 
the transmission matrix of the structure being considered 
with the scattering matrix approach, which is relatively 
easy to implement and exhibits good numerical stability, 
also in the case of high values of the magnetic field B. 

In particular, we consider a 2-dimensional Schrodinger 
equation in the x-y plane, with x being the direction 
of electron propagation (and thus z the direction of the 
uniform orthogonal magnetic field). We choose a vec- 
tor potential with a single nonzero component along the 
longitudinal direction x (with a value —By) and we sub- 
divide the structure into a number of transverse slices, in 
each of which the scalar (confinement) potential can be 
assumed to be independent of x. In each slice we expand 
the transverse eigenfunctions onto a basis made up of the 
transverse eigenfunctions for B = 0; the coefficients of 
this expansion and the associated longitudinal wave vec- 
tors are found solving an eigenvalue problem llfill . Then, 



with a mode-matching technique, we compute the scat- 
tering matrices of the sections extending from the middle 
of each slice to the middle of the following slice (and thus 
containing only one discontinuity of the potential) . Com- 
posing the scattering matrices of all the sections m , we 
find the S-matrix (and, as a submatrix, the transmission 
matrix t with elements t nm ) of the overall structure. 

The value of the conductance and of the shot noise 
power spectral density can be then computed by means 
of the relations 0 

g = y; \t nm \ = —j^- (i) 

n,m j 

and 

S i =A ~h ( 2 ) 
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where the Wj’s are the eigenvalues of the matrix tt' , e 
is the electron charge, h is Planck’s constant, and V is 
the externally applied voltage. Since the power spectral 
density of full shot noise is given by Sj s = 2e\I\ (I being 
the average current through the device) the Fano factor 
7 can be written as 



7 = 



,r /) 



(3) 



For a small number of modes propagating through the 
constrictions, quantum interference effects lead to wide 
relative fluctuations, as a function of energy, of both the 
numerator and the denominator of this expression. Av- 
eraging is therefore needed, and care must be taken to 
perform it correctly: in order to measure shot noise, the 
applied voltage V must be much larger than kT /e ( k be- 
ing the Boltzmann constant and T the absolute tempera- 
ture) . Based on the detailed expression of the shot noise 
power spectral density at finite temperature provided by 
Biittiker E3, in the case of eV kT a good approxima- 
tion of the shot noise term is given by a uniform average 
over an interval eV around the Fermi level of eq. An 
analogous conclusion, i.e., uniform averaging over the in- 
terval eV, can be drawn for the conductance term. We 
remark that it is essential that averaging be performed 
separately on the numerator and the denominator before 
taking the ratio (as in the actual measurement proce- 
dure). 

Each of our data points has been obtained by averag- 
ing over the noise power spectral density and the con- 
ductance for 41 different energy values in an interval of 
width 0.22 meV around the Fermi energy. In all of the 
following calculations we consider a Fermi energy Ef of 
9.134 meV and the effective mass of gallium arsenide, 
m* = 0.067mo, with mo being the free electron mass 
(values consistent with the experiments by Oberholzer et 
al.). 

This numerical approach has been used for the simula- 
tion of shot noise suppression in a rectangular cavity, first 
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FIG. 2: (a) Fano factor as a function of the magnetic field 
computed for a mesoscopic cavity with constrictions of 40 nm 
(solid circles), 60 nm (empty circles), 100 nm (solid squares), 
300 nm (empty squares). The diamonds represent the experi- 
mental data by Oberholzer et al; (b) same data as in (a), but 
plotted as a function of the ratio of the classical cyclotron 
diameter Dc to the constriction width Wc- 



in the absence of a magnetic field. In particular, we show 
that the expected shot noise suppression is achieved with 
a simple hard-wall model defining a rectangular cavity, 
5 //in long and 8 //m wide (corresponding to the litho- 
graphic dimensions of the cavity in refs. 00) , with 
symmetric constrictions whose width is chosen on the 
basis of the desired number of propagating modes. 

As constrictions are made wider (reaching a condition 
in which the analytical theory yielding the value 1/4 for 
the suppression factor is not applicable any longer), we 
observe a decrease of the Fano factor, until it drops to 
zero as their width equals the cavity width and the struc- 
ture becomes a noiseless perfect quantum wire. In the 
lower panel of fig. Q we report the Fano factor from our 
numerical calculations (empty circles) as a function of 
the number of propagating modes: empty squares rep- 
resent the experimental data obtained by Oberholzer et 
al. 0 , and the agreement appears to be good, without 
the need for any fitting parameter. In the upper panel of 
fig. ID we report the ratio of the overall cavity resistance 
to the sum of the constriction resistances, as a function 
of the number of propagating modes with a solid line 
for our rectangular cavity and with a dotted line for a 
stadium-shaped cavity (8 //nr wide and 9 ft m long, since a 
horizontal stadium cannot be made with the same aspect 
ratio as that of the experimental cavity): such a quantity 
drops below unity, in disagreement with the experiment. 
We show also the calculation for the stadium-shaped cav- 
ity in order to exclude the possibility of a significant role 
played by a classically chaotic shape. Also the inclusion 
of disorder in the cavity (data for which are not reported 
here, due to space constraints) would not provide a valid 
justification of resistance behavior, since it would raise 





FIG. 3: (a) Fano factor as a function of the cavity width 
W (solid circles) and length L (remaining plots) for 100 nm 
symmetric constrictions; (b) conductance through the cavity, 
normalized with respect to the conductance quantum, as a 
function of the orthogonal magnetic field, for two values of 
the constriction widths. 



both resistance and noise. 

We suggest, as a possible explanation, that such an 
additivity over the whole measurement range can be the 
result of partial electron thermalization: as constrictions 
are made wider and electron diffraction is reduced, the 
total resistance of the cavity should start decreasing be- 
low the sum of the individual constrictions, according to 
the outcome of our ballistic calculation, while it would 
remain constant, independent of constriction width, if 
the cavity were a real reservoir, with electrons in ther- 
mal equilibrium with the lattice. If full thermalization of 
the electrons did occur, shot noise would also disappear 
(as long as the constrictions have integer transmission); 
a partial thermalization would instead determine an in- 
crease of the resistance toward the classical sum rule and, 
at the same time, a decrease of the noise, with respect to 
the ballistic result: this is consistent with the comparison 
between our results and those from the experiment. 

We have then moved on to the simulation of the be- 
havior of a cavity in the presence of a magnetic field 
orthogonal to the plane of the device, with a value up to 
a few tesla. 

Results for the Fano factor are reported in fig. Eta) 



4 



for symmetric constrictions with a width of 40 nm (solid 
circles), 60 nm (empty circles), 100 nm (solid squares), 
and 300 nm (empty squares). The Fano factor exhibits a 
decrease as a function of the magnetic field, with a rate 
depending on the width of the constrictions and, there- 
fore, on the number of modes propagating through such 
constrictions. The smaller the number of propagating 
modes the larger the fluctuations, due to the more sig- 
nificant relative contribution from quantum interference: 
for the narrowest constrictions residual fluctuations sur- 
vive our averaging procedure. 

The only available experimental data, from ref. 0, 
are for a single mode propagating through the constric- 
tions 20], which corresponds, in our model, to a con- 
striction width of 40 nm. Experimental results from 40] , 
represented with diamonds, lie slightly below the theo- 
retical data for the 40 nm constrictions and close to those 
for the 60 nm constrictions. Thus our simulations tend to 
slightly overestimate the Fano factor, for narrow constric- 
tions and large B values. This can be attributed, in par- 
ticular for a single propagating mode, to the specific scat- 
tering properties of hard-wall constrictions, which tend 
to create a set of modes inside the cavity, some of which 
are almost completely reflected at the exit constriction, 
contrary to gradual adiabatic constrictions (as we have 
verified with numerical simulations on a structure with 
openings defined as in ref. pljl. 

In fig. [lb) we report the same data as that of fig. Ha), 
plotted as a function of the ratio of the classical cyclotron 
diameter Dc (twice the cyclotron radius) to the constric- 
tion width Wc'- it is apparent that results for different 
constriction widths are essentially superimposed, thus 
demonstrating our conjecture that Dc/Wc is the actu- 
ally relevant parameter defining the level of shot noise 
suppression, simply because it determines the amount of 
scattering, and therefore diffraction, that occurs at the 
constrictions. Indeed, also the r.h.s. of eq. (4) of ref. 0] 
is proportional to Dq and inversely proportional to Wc'- 
however it would also be proportional to the perimeter 
(or length), a dependence which is not observed in our 
results. 

Running our simulations for different values of the cav- 
ity width and length (see fig. Ha)), we have noticed that, 
in the regime of narrow constrictions and noise suppres- 
sion due to the magnetic held, the Fano factor has only 
a mild and nonmonotonic dependence on such parame- 
ters (unless the constrictions are rather wide, and we are 
not in the regime of interest any longer), arguably re- 
sulting just from fluctuations due to interference effects. 
This is not in agreement with the interpretation of the 
phenomenon proposed by Oberholzer et al. 0 . Start- 
ing from a first-order expression involving the ratio of a 
quantum scattering time to the dwell time in the cavity, 
they suggest that the Fano factor should depend on the 
cavity area available for electron motion and they explain 
the shot noise suppression with a reduction of this area 
as a consequence of the formation of cyclotron orbits. 

From the results of our numerical calculations, we 



propose a different and intuitive interpretation, directly 
based on the comparison of the cyclotron diameter with 
constriction width. In a semiclassical picture with skip- 
ping orbits crawling along the cavity walls, if the cy- 
clotron diameter Dc is much larger than the constriction 
width Wc, an electron impinging against such a constric- 
tion is likely to be reflected, thereby leading to diffrac- 
tion. On the contrary, if Dc is smaller than Wc, it is 
very likely that the electron traverses the two constric- 
tions without undergoing reflections, so that shot noise 
is strongly suppressed (in the limit of no scattering, the 
Fano factor would drop to zero). 

We have also computed the conductance of the cavity 
as a function of magnetic field, observing the transition, 
measured in ref. (l3j | . between a condition, for B = 0, 
in which it equals one half of the conductance of each 
constriction, and that for a large magnetic held, in which 
edge states are formed and it reaches the value of a single 
constriction. Results are reported in fig. Hb), where the 
conductance for a constriction width of 40 nm (1 prop- 
agating mode) starts at 0.5 units and increases toward 
1 unit as in the experimental results of ref. 10, while 
the curve for Wc = 100 nm (4 propagating modes for 
B = 0) has an initial value of 2 and a value of 3 for a 
large magnetic field. This latter value is not 4, as ex- 
pected, because the total number of modes propagating 
through the constrictions decreases down to 3 at about 
0.5 T. We notice that the experimental data for conduc- 
tance are in this case more closely reproduced by our 
ballistic simulation, as can be explained by the fact that 
in the presence of a magnetic field the effective area of 
the cavity and the time spent by each electron inside it 
are reduced, and therefore thermalization is decreased. 
We also wish to point out that, for the explanation of 
the observed phenomena, there is no need to invoke the 
presence of irregularities or defects in the cavity: indeed, 
as in the case of no magnetic field, they would lead to an 
increase of noise, due to the resulting electron diffraction. 

In conclusion, we have presented a numerical simu- 
lation of shot noise suppression in mesoscopic cavities, 
reproducing the results of recent experiments without 
the need for fitting parameters. In particular, we get 
very good agreement with the measured behavior of the 
Fano factor as a function of the width of the constric- 
tions and, although with fluctuations due to interference 
effects, we are able to closely reproduce the dependence 
of the Fano factor on magnetic field. From our results, 
we observe that the reduction of the Fano factor with 
increasing magnetic held can be simply explained as the 
consequence of the decrease in the ratio of the cyclotron 
diameter to the constriction width, which leads to pro- 
gressively improved transmission through the apertures, 
and therefore suppression of diffraction. 

This work has been supported by the Italian Min- 
istry for Education, University and Research, through 
the FIRB project “Nanotechnologies and nanodevices for 
the information society” and through the PRIN “Excess 
noise in nanoscale devices.” 
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